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o 

^vj , Abstract. In this paper we study the heat trace expansion of the perturbed 

■ harmonic oscillator by adapting to the semiclassical setting techniques devel- 

J"~H' oped by Hitrick-Polterovich in [HP]. We use the expansion to obtain certain 

~~ ' inverse spectral results. 
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PLh ! 1- Introduction 

C/3 ' 

Hitrik and Poltcrovich obtained in [HP] a simple formula for the on-diagonal 

heat kernel expansion of the Schrodinger operator, —A + V, with V £ C°°(K™) a 

bounded real- valued potential. In this paper we apply their techniques to study the 

semiclassical behavior of the on-diagonal heat kernel expansion for the perturbed 

semi-classical harmonic oscillator 

^\ . More precisely, we consider the kernel of the operator e but with ordinary time 

OQ ' t replaced by 



1 , 1 



2 : (L2) f ^ log i-^ = 2 H 1 



which greatly simplifies the calculations in the semiclassical regime. 
Our first result is: 



Theorem 1.1. Assume that V £ C°°(IR' 1 ) is bounded below and that it and all its 
derivatives have at most polynomial growth at infinity. Then, on the diagonal, the 
Schwartz kernel of the operator e~ tH , where t is given by (1.2(1 . has an asymptotic 
expansion as H tends to zero of the form 

(1.3) ft 2 5> 2fc T fe ( s ,x). 

fc=0 

Moreover, the first three coefficients in this expansion, integrated over W l , deter- 
mine the following guantities: 
(1.4) 

V(x)e- S ^ 2 dx, [ V 2 {x)e~ s \ x \ 2 dx, [ (V 3 (x) - VAV) e^ 8 ^ dx. 
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The Hitrik-Polterovich method, adapted to the present situation, results in a 
procedure to compute the T&. In particular To(s, x) = 2sV(x)e~' x ' s . 

For ft sufficiently small the spectrum of H is discrete and the quantities (|1.4|l arc 
spectral invariants of V associated to the ^-dependent spectrum of H. By analyzing 
these invariants we obtain several inverse spectral results, namely: 

Corollary 1.2. Let S r = {a; 6 M" ; |x| = r}. The following properties of V can be 
detected from the h-dependent spectrum of H : 

(a) Whether V is constant on a given sphere S r , and if so the value of the 
constant. 

(b) Whether V is compactly supported, and if so the the smallest annulus about 
the origin containing the support of V . 

(c) Within the class of odd functions V, one can determine whether the restric- 
tion of V to any annulus about the origin is linear on that annulus. 

Remark 1.3. Item (c) is a consequence of a much stronger but slightly more tech- 
nical result (see Proposition EH ) 



The basic ingredients in the derivation of the expansion (JT73J) are a variant of 
Mehler's formula and the Kantorovitz formula for expressing the heat expansion of 
the sum of two operators, A and B, in terms of the heat expansion of B alone. (The 
latter is also the basic ingredient in the proof of the Hitrik-Polterovich result.) We 
will discuss Mehler's formula in §2 and the Hitrik-Polterovich formula in §3. Then 
in §4 we will describe what this formula looks like if one replaces the standard heat 

\ X _ y I 2 

kernel, (47ri)~2"e ** , by the semi-classical Mchlcr kernel. As mentioned, the 
expansion of this formula in powers of h 2 generates a sequence of heat trace invari- 
ants, and in §5 we will discuss a symbolic method for computing these invariants. In 
§6 we will illustrate these methods by computing the first three of these invariants, 
and finally in §7 we will prove the aforementioned inverse spectral results. 

2. Mehler's formula 
Let L be the operator 

if- 



^ L =2"£(-^ 



,.2 



Mehler's formula for the Schwartz kernel of e tL is 

(2.2) 

e~ tL (x,y) = 7r~^(l-e -2 *) -: J exp <^ — - 



1-e 



-2/ 



\v\ 2 n ,~2t 



(1 + e~ M ) - 2e~ t x ■ y 



(see for instance [Simj page 38). Rescaling the variables x and y by the factor 1/vH 

and t by the factor ft we get, for the heat kernel of the semi-classical harmonic 

oscillator 

ft 2 A \x\ 2 nh 

A = A + — , 

2 2 2 ' 
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the expression 

(2.3) 

e-"(s, y) = .-f (l-e-^)-f exp |-__L- H - 

The term in square brackets can be rewritten as 



\^l±^l(l + e- 2ht )-2e- th x-y 



^-(l-e- tH ) 2 + e- m \x-y\ 2 , 



2 
and hence the term in curly braces is equal to 

(24 ) M 2 + | y | 2 l- e -«» ^ 

1 ' J 2ft l + e- ts ft(l-e- ift )(l + e-* a )' *" 

Now introduce the new time scale 

(2-5) S = ^TtS = ^ (1 + 0(<2 " 2)) 

or, alternatively, 

(2.6) e -^=i^f. 

v 7 l + hs 

Then the expression (|2.4[) becomes 



x\ 2 + \y\ 2 if 11 , a 



2 4U 2 s 

or 

/ 7 n k-y| 2 N + yl 2 

and hence for the Schwartz kernel of heat operator e~ tA we get the formula 
(2.8) e~ tA {x, y) = (47rfi)" * a"* (1 + sft) n cxp f-^r^ 

3. The Kantorovitz formula 

Let A and £? be linear operators on an appropriately defined Hilbert (Banach, 
Frechct, • • ■ ) space which generate strongly continuous semigroups e tA and e tB and 
such that the sets of C°° vectors satisfy: D°°(A) dD oa (A + B). (Both conditions 
will be satisfied automatically in what follows.) Then according to Kantorovitz, 
jKaj . e^ A+B ' can be expressed as a series 

(3.1) e*( A+s > = (/ + tX x + jX 2 + ■ ■ ■ )e tA , 

where the Xi's are defined by 

(3.2) X =1, X 1= B, X 2 = B 2 + [A, B] 
and in general 

(3.3) X m = BX m -i + [A, X TO _iJ. 
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There is also a simple closed form expression for X m : Letting H = A + B, 
(3.4) X m = H m - mH m - 1 A + ( m ) H m ~ 2 A 2 + ■■■ . 



2, 
Example: ([HP]) Let A = -A R „ and H = -A R „ + V. Then 

(3.5) e- tA (x,y) = (4irt)-%e- lE ^ jl 
and hence 

(3.6) e~ tH (x, y) = (47rt)"9 V (-l) m — f a m (x, y, t), 



m=0 



where 



(3.7) a m (x,y,t)= £ (^Va^-^. 

l+j=m ^ ' 

From this identity it is relatively easy to get an asymptotic expansion of e~ tH (x, x) 
as a Taylor series in t for which the summands are (at least in principle) computable. 
(See [HP], §2 for details.) 

In the spirit of this example, let 

„ ^ / ft 2 d 2 x 2 ft 

I— 1 v * 

and let H = A + ft 2 V. Then, as above, 

(3.9) e~ tH (x, y) = V (-1)"*^ o^ar, y, i, ft), 

* — ' to! 

m=0 

where 

(3.10) a m (x,y,t,h)= £ (-1)^ (^F^e^^y). 



By (|2.8j) the computation of this sum reduces to computing 
(3.11) £ (-l^W'exp (-^f ) /(*,»,-)«' 

where s is given by (|2.5[) and 



(3.12) f(x,y,s) = expl -a 

The expression above is similar to the Hitrik-Polterovich expression 

i+ E (7)(-A + F )i A^ P (-Ji^!!)r, 

except for the presence of the factor (|3.12p . However, since we'll mainly be inter- 
ested in the ft dependence of the expression p. lip , and p,12j) depends in an explicit 
way on ft, our computations will be very similar to theirs. 
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4. Computations 
As above let A be the operator 



El h 2 d 2 x% n \ 



Jdx 2 ' 2 2 

B the operator, h 2 V, and X m , m = 0, 1, 2, • • • the operators defined inductively by 

(4.1) l m = BX m _ 1 + [A,l m _ 1 ] 

and Xq = I. It will be convenient to write this formula as 



(4.2) 



A x 2 

X m = fi, VX m -i — h [— ,X m -i] + [— ,X m -i] 



From this formula one gets: 

Proposition 4.1. The operators X m are of the form 



(4.3) 



x m — h m 2__, h l Xm , 



i=l 
i=m mod(2) 



where X l m 1 is a differential operator of degree i — 1 not depending on h. Moreover, 
these operators satisfy 



(4.4) 



X 



m+l 



i-1 



-i X m 



-,x 



i+1 



vx. 



(-1 



The proof is a simple inductive argument. 

To compute the m th summand in the Kantorovitz expansion (|3.ip . we must 
apply X m to the Mehler kernel 



(4.5) 

and then set x = y 
Proposition 4.2. 



\%-y\ 2 \x + y\ 2 

e(x, y, s, K) = exp ( -^- s 



X m (e(x,y,s,h))\ x=y = h m Y^ ^Ki l { e {x iy ,s,h))\ x=y 



i=l 
i=m mod{2) 



is equal to: For m odd and with I = m 2 , 



h m+1 \Y^e m , r {x,s)h 2r \ 



s-'e— I*' 



(4.6) 

and for m even and with 1=^ — 1, 

(4.7) 

where in all cases the e m ^ r are polynomials in s of degree at most 2r. 



h m+2 (j2e„ hr (x,s)h 2 A 



a -i e ->\*\' 
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Proof. We first note that for multi-indices, /z, 

(4.8) <9£e~^S 



x=y 



Cjl m S 2 



for even p and for non-cven p, where c M = (— t)' 1 ''^ for p = 2v. The result 
follows from this, Leibniz' formula, and the properties of the operators X£^ • 



□ 



Thus making the substitution 

, , l , i + hs ( h 2 s 2 h 4 s 4 

(4.9) t = - log — — - = 2s 1 



ft 1 — fts \ 3 5 

the t m X m e(a;, y,i, ft)|a; = y term in the Kantorovitz formula gets converted into 

(4.10) 2 m h m+1 s l+1 [J2e m Ax,s)h 2r ) (l 



\r=0 



I = m j^ 1 , for m odd, and 



(4.11) 2 m h m+2 s l+2 J2 e "^( x > «)ft 2r ( 1 



\r=0 



1=^ — 1, for m even. 



5. Symbolic features of the expansions (|4.6|) - (|4.7|) 
We showed above that there exist functions p m (x, s), m = 0, 1, . . . such that 

(5.1) X m e(x,y,s,h)\ x=y = H m+1 p m (x,s) + 0(h m+3 ) 
for m odd and 

(5.2) X m e{x,y,s,h)\ x = y = H m+2 p m (x, s) + 0(ft m+4 ) 

for m even. We will show in this section that for m odd, p m (x,s) is computable 
purely by "symbolic" techniques and will prove a somewhat weaker form of this 
assertion for m even. 
Let 

(5-3) X l m+l = J2 <rn+l(x)D a 

\a\<i 

and let 

(5-4) ?4 + l = E <m + l(x)C 

be the full symbol of X % m+1 From (|4.4[) and standard composition formula for left 
Kohn-Nirenbcrg symbols one gets 

(5.5) 

*U - f E ^f ^- - ^J ^+ (E ^^ + ^ J p^+n*^ 1 - 
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In particular, if 

<J m+l = 2-^i a i,m+l\ x fc 

\a\=i 

is the principal symbol of X 1 nl+1 and 

a m+l = 2-^/ a i,m+l\ X )€ 

\a\— i— 1 

the subprincipal symbol, we get from (|5.5[) that 



(5.6) 

and 

(5.7) 

~i 

Letting 



'm+l 



1 y> ( t d i_ x _d_ i+1 \ 



r—1 v / r—1 x ' 



Wr ^\ + v^ 



<? m = h m £ ^v;- 1 , 



z=m mod(2) 



and letting IA be the raising operator 

WftV = h l+2 a 
for i > 0, we can write these formulas more succinctly in the form 



(5.8) 

and 
(5.9) 



^gfc^-^r" 






13 2 „ 1_92 
29^ 



CTm-l+VKa 



m-li- V UUm-1- 



d .. _d_ 

Xr ~Ew r 
d£, r 



h 2 v 



In particular, iterating (|5.8|) we get 

(5.io) „ m = E ^^/. 

v r—1 v 

and, as special cases of (J5.1QD , 

(5.1D T 1 - fattr&y 1 * 

As applications of these formulas let m be odd and consider the zth summand of 



X m e(x, y, s, h)\ x=y = h m Y^ ffXln ^(x, y, s, h) 



i=l 
i=m mod(2) 
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-y\- 



\a\= i— 1 



— s\x\ 



x=y 



By ()5.3j) this is equal to 
(5.12) 

plus terms of order 0(h m+3 ) and by (|4.8j) . (|5.12[) is equal to 

(5.13) h m+1 i J2 Ci, m wJs-¥ e -N ! 

To summarize: 

Proposition 5.1. For to odd, the leading term p m in 115.1]) is given by 

rn 

Pm (x,s) = e- S ^ 2 ^ S ^ Yl a i-hm( X ) C c- 

i=l \a\=i— 1 

i=m mod(2) 

For each i, the quantity Y]i a |— f— i. a f-i ml 3 ') " * s obtained from the principal symbol 
of X^ 1 by substituting every monomial £ Q by the constant c a . 

For to even the computation above is similar, however one gets h m+2 contribu- 
tions to (15.21) from both the terms 



h m + i J2 a%_ hm {x)D<*e 

\a\=i-l 



\*+vV 



x=y 



and the terms 



hm+l J2 aU, m {x)D a e-^ 

\a\=i-2 

The second summand (involving the subprincipal symbol of -X^i ) is as before, 
(5.14) h m + 2 £ CLJ^^^e-^l 2 , 

\a\=i-2 

but the first summand (involving the principal symbol of X^ 1 ) becomes 



(5.15) 



nm+ 2 J2 a?_i, m (x)c al r )a - 



2 x r e 



-s\x\ 2 



\a\— i — 1 
l<r<n 



where a^ r > — (ai, • ■ • , a r — 1, • • • , a„). This 



proves: 



Proposition 5.2. For to even t/ie leading term of li5.S\) depends only on the prin- 
cipal and subprincipal symbols of X m . 

We now explore some spectral consequences of the previous results. 

Proposition 5.3. For to odd the quantities J p m (x, s) dx are spectral invariants of 

V. 
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x a irw e S ' x ' ' w here \a\ + \j3\ < m — 1. The associated contribution to the heat 



Proof. For to odd we can, by (|5.10[) . express p m {x, s) as a sum of terms of the form 

trace 

f d^V 

< 5 ' 16 > J *°W> 

can, by integration by parts, be written as sums of integrals of the form 

jW-i^, M < m -i. 

Thus 

(5.17) / p m (x,s)e- s]xl2 dx= I ' P {x,s)Ve- s ^ 'dx, 

where p{x, s) is a universal polynomial of degree m — 1 in x. Moreover, for every 
A € SO(n) the heat trace expansion for the potentials V and V A , where V A (x) = 
V(Ax), are the same. Hence by averaging over SO(n) we can assume that p(x, s) 
is SO(n) invariant, i.e. 

fc TO- 1 

p{x,s) = ^2 Xl {s)\x\ 2 \ k = ^~ 

8=0 

and thus (|5.17[) becomes 

//•oo k „ 

Pm (x,s)e- S ^ 2 dx= drY,Xi(s)r 2l e- sr l V(x)da r , 

JO i=1 J\x\=r 

where da r is the standard volume form on the (n — l)-sphcrc \x\ = r. Wc will sec 
below however that, for each r > 0, the integral 



(5.19) / V(x)da r 

J |x|= r 

is itself a spectral invariant of the perturbed harmonic oscillator and hence the 
terms K' n+1 J p m {x,s)dx in the heat trace expansion above can be read off from 

it. □ 



In the case m even one also gets a similar description of the contributions of ([57 
to the heat trace. The contribution coming from the term (|5. 15[) only depends on 
cr^" 1 and hence as above is expressible in terms of (|5.18|) : and as for the contributions 
coming from (|5.14[) one can prove by induction that these give rise to heat trace 
invariants which are sums of expressions of the form (|5.16|) and of the form 

(5.20) f x a V^er s ^ 2 dx. 

J ax p 

Indeed the second summand in (|5.7p is purely symbolic; so as we've just seen it 
contributes terms of type (|5. 16[) to the heat trace. Similarly the third summand 
contributes terms of type (|5.20p and by a simple induction on m one can show that 
the first summand of (|5.7|) is a linear combination of terms of the form 

dx~< dx 5 
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and 



Cx p 



dW 

dx~< ' 



These give rise to contributions to the heat trace of the form (|5.16[) and 

u — — e s]x] dx, 

ox v ox 1 

which by integration by parts can be written as expressions of the form (|5.20l) . 
Finally, the 0(n) invariancc of the heat trace enables one to simplify these further 
and rewrite them as sums of the form 

J^V^x^T^Vdx. 

6. The first heat invariants 

It is easy to see (either by direct computation or by the symbolic formulas in the 
preceding section) that for m < 4 the X^ s are given by 

(6.1) X° = V, 



(6.2) 








(6.3) 

Y 2 


E 


d 2 V 


d 2 


A 3 


dxidxj 


dxidxj 



xl 



v^ dV d T o AV 
*-" dx, dx, 2 ' 



yl ( A7-V)A + ^-V. 

Z^ Fln-.- K ■) > J Fin-.- A 



dXi 



dx, 



■V J 



VAV 



(6.4) 
(6.5) 



A 3 ~ 2^ Xl Fln-' 



' dxi 



X 



A ^ d 2 V 



if- 



2 ' ^-^ dxidxj dxidxj 

plus terms of degree less than two, and 
(6.6) 



_d_ 

dxi 



££-<^-^ 2 >£-i+I> 



0_ 
dx,' 



2 V d 2 
dxidxj dxidxj 



X 



•E< 



dV 
dxi 



x"_ ^ d 2 V 



d 2 



d 



dxidxj dxidxj 



£^ Ay 



2 ' dx, 



E 



xt-g-V*. 

dxi 



Thus the h 2 term in the heat trace expansion determines 
(6.7) fve- s ^ 2 dx, 

and hence by the inverse Laplace transform determines the integral 

Vdcr r . 



for each r > 0. 

The h A term involves the m = 3 contribution of (|5.1|) . but as we saw above this 
is expressible in terms of (|6.8|) . As for the contribution of (|5.2j) to the fi 4 term, the 
first and third summands can be converted by integration by parts into integrals 
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which are expressible in terms of ()6.8[) and the second summand gives a new heat 
invariant, 

(6.9, /^e-NlV, 

which, by the inverse Laplace transform, is convertible into 
(6.10) / V 2 da r . 



The h term in the heat trace expansion involves the m = 5 contribution of (|5.1[) 
which, as we saw in the previous section, is expressible in terms of (|6.8p . the m = A 
contribution of (|5.2[) . which is subprincipal and hence only involves the terms in 
(|6.6|) (all of which can be converted, by integration by parts, into expressions in 
(|6.8[) and (|6.10[1 ). and the cubic and quadratic terms in (|6.5|) all of which, except 
for the term 

^ d 2 V d d 
^—^ dxidxj dxi dxj ' 

can be converted by integration by parts into expressions in (|6.8[) and (|6.10p . Finally 
the H 6 terms coming from (|6.3[) and (|6.4[) are all convertible by integration by parts 
into expressions in (|6.8I) and (|6.10[) except for the last summand of (J6T3J): the term 

(6.12) V-™ 

The term (|6. 1 1[) gives, by (j5.14[) and (|4.11[) a contribution 

(6.13) _ ft V^e-W a 

to the heat trace expansion, and the term (|6.12[) gives, by (|4.10p . a contribution 

(6.14) h e s 2 (V 3 - -^) e - s l-l 2 

to the heat trace expansion; hence the sum of these two terms gives rise to a new 
heat trace invariant 

(6.15) f(V 3 ~VAV)e- sM2 dx 

which, by the inverse Laplace transform, can be converted into the invariant 

(6.16) l {V 3 -VAV)da r . 

J \x\—r 



This finishes the proof of Theorem 11.11 

7. Applications to inverse spectral problems 

In this section we apply the first heat invariants (|6.8j) . (|6.10|) and (|6.16|) above 
to the inverse spectral problems of recovering information about V from the H- 
dependent spectrum of H. 
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Fix any r > 0. First let's consider V which minimize the second invariant (|6.10|) 
subject to the constraint 

(7.1) / Vda r = constant. 

J a:|— r 

According to the Cauchy- Schwartz inequality, the minimizers are exactly those 
functions V that are constant on the sphere \x\ = r. It follows that the set of 
potentials V that are constant on a given sphere \x\ = r is intrinsically defined by 
its spectral properties. Moreover, one can spectrally determine the constant value 
for each potential in this set. This proves parts (a) and (b) of Corollary 1 1.2 1 

Next let's assume that V is an odd potential. (Using band invariant techniques, 
one can show that being odd is also a spectral property, c.f. |GUWj .) As we have 
seen that the invariants determines ||V"||^ 2 (g ■> for each r. Taking the r derivative 
of the second invariant, we get that (V, ^r)z, 2 (5V) is a spectral invariant. On the 
other hand, since V is odd, the third invariant (|6.15[) becomes 



(7.2) - / VAVe 



Recall that in spherical coordinates 



AT , d 2 V n-ldV 1 . „ 

AV = ^TT + a- + ^ A SrV- 

or z r or r^ 

A simple computation shows that 

v n^ldV_ sr%n _ ldrd ^ = n-1 r 2 d L- sr%n ^\ drdGr 
r or 2 J dr V / 

which is a spectrally determined quantity, since we know the integrals (|6.10[) for all 
r. Similarly, 

V%X-er sr \ n ~ X drda r = f (^-) e^ 'r n - 1 drda r + I ' V^-^-( e - sr ' 'r^drdar, 
or z J \ or J J or dr 

and again the second term on the right 

V^-^{e- sr \ n - x )drda r 
or dr 

is also spectrally determined according to (|6.10p . It follows from (|7.2| that the 
integral 

(7.3) I yjf) e ' sr2rn ~ ldrd(J r- f \vA Sr Ve- 8r2 r n - l drd(x r 

is spectrally determined, which, by the inverse Laplace transform, gets converted 
to the invariant H^-H^g \ + ^{V,~A Sr V) L 2( Sr ^ for each r > 0. Note next that 
for every V one has the following inequality 



\nl H s T) {\\%\\l> { s r) + l -(V,-A L , iSr) V) L . {Sr ^ > (V,^)i HSr) +±\\V\\ 



£ 2 (Sr)' 
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where Ai is the first eigenvalue of the (non- negative) Laplacian on S r . Both sides 
of the inequality are spectral invariants, and equality holds if and only if V satisfies 
the conditions 

dV 

(7.4) -rr~\s = X v \s an d V\s is a spherical harmonic of degree one 

or 

where x is a constant, and if so one can determine X- This proves the first part of 

the following 

Proposition 7.1. The class of functions defined by the conditions |7.^| ] is spectrally 
determined. Moreover, for any V in this class one can determine the ratio 

dV IV 

on a given sphere S r . 

The determination of \ (which of course can depend on r) is done by looking at 
the quotient of (V, ^>L»(s r )/l|V||| a(Sp) . 

As a consequence, one can determine whether a potential is of the form 

V(x) = f(r)g(cr) 

on a given annulus r\ < \r\ < r%, where g(o~) is a spherical harmonic of degree one, 
and if so, determine the function f(r). In particular, one can spectrally determine 
linear potentials on any annular region r\ < \r\ < r 2 : They are just the potentials 
in the previous class with x = r - 
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